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ABSTRACT 


The circulation in the South China Sea (SCS) is mainly determined by the monsoons. 
This monsoon-induced circulation is connected with the circulation in the Pacific Ocean to 
a small extent and linked with that of the Indian Ocean only by the Timor Current. During 
winter, the water enters the South China Sea at the north through Formosa Strait and Luzon 
Strait, and exits at the south through Karimata Strait. During summer, the circulation re- 
verses, the water enters at the south and exits at the north. Water from the north 1s generally 
cooler than water from the south. The northward (southward) flow induces warm (cold) ad- 
vection. Therefore, the seasonal cycle of the circulation in the south China Sea becomes a 
major factor controlling the seasonal cycle of the sea surface temperature, which might ef- 
fect the seasonal variation of the Eastern Asian monsoon system. 
In the preliminary results, the seasonal circulation and temperature patterns simulated 
by the primitive equation model are similar to those reported by Wyrtki (1961). Yet, our 


model 1s better than the Pohlmann (1987) shallow-water model and shows more details. 
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I. INTRODUCTION 


The South China Sea, with the Gulf of Thailand included, is the largest marginal sea 
of the Southeast Asian Waters (Figure 1). It is separated from adjacent oceans by the chain 
of the Southeast Asian Islands and contains the broad shallow Sunda Shelf on the south, a 
shelf extending from the Gulf of Tonkin to Formosa Strait on the north, a deep sea basin in 
the center, and the Dangerous Ground with numerous reef islands over the southeast 
(Figure 2). From 1959 to 1961 a Naga expedition was carried out by the Sripps Institution 
of Oceanography of University of California and jointly sponsored by the governments of 
the former South Vietnam, Thailand, and the United States of America. It provided the 
observational results of marine investigations of the South China Sea (Wyrtki, 1961). 

In summer (June), the southwest monsoon prevails and surface currents flow 
northeastward entering the South China Sea through the Gasper and Karimata Straits 
(Figure 3A). The largest part of the South China Sea water masses pass south of Taiwan 
into the root of the Kuroshio Current; a small part flows northwards through the Formosa 
Strait joining the Kuroshio Current. A current develops off the coast of Borneo towards the 
southwest with a current speed between 0.12 m/s and 0.25 m/s, which turns back to the 
north in the region of Natuna Island. Along the coast of Luzon, water flows northwards 
from the Sulu Sea with a surface current speed between 0.12 m/s and 0.25 m/s. A jet current 
flows from Vietnam to the Luzon Strait with a current speed of 0.25 m/s to 0.30 m/s. 

In winter (December), the northeast monsoon prevails and surface currents flow 
southwestward entering the South China Sea through the Formosa and Luzon Straits 
(Figure 3B). Waters from the Kuroshio Current and the East China Sea enter the South 
China Sea flowing southwestward along the coastline of Asia and flow out of the South 


China Sea through the Gasper and Karimata Straits. The surface current speeds are about 
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Figure | : Location of the South China Sea 
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Figure 3A ; Observational Summer (June) Surface Circulation of the Southeast 


Asian Waters (Wyrtki, 1961) 
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Figure 3B : Observational Winter (December) Surface Circulation of the Southeast 


Asian Waters (Wyrtki, 1961) 


0.5 m/s off the Hainan coast, about 0.75 m/s off the coast of Vietnam, and about 0.3 m/s in 
the Gasper Strait. The counter current in the center of the South China Sea is present and 
lies farther to the east because of the strength of the main current. From the Sulu Sea water 
flows into the South China Sea and is mixed into the main current. A current boundary 
exists between Vietnam and Borneo and a very weak (less than 0.12 m/s) southwestward 
flow occurs along the coast of Borneo which encounter the waters from the north flowing 
northward along the right hand side of the current boundary. 

Pohlmann (1987) used a 12-layer shallow water circulation model to study the 
responses of the South China Sea to the east Asian monsoons. The model horizontal grid 
size is about 50 km. The simulated mean surface circulation for the summer and winter 
seasons is shown in Figures 4A and 4B, respectively. Since they do not show the current 
speed clearly, we only describe the surface circulation patterns as follows. 

In summer, waters from the Java Sea enter the South China Sea through the Gasper 
and Karimata Straits flowing clockwise and eastward from the southern coast of Vietnam 
to Borneo. Then, it flows northward along the coast of Palawan and Luzon to west of the 
Luzon Strait. In the northern South China Sea, there is a cyclonic circulation between 
Hainan and the Luzon Strait. A small amount of Kuroshio current flows into the South 
China Sea and meets the waters from Borneo flowing westward, and then it encounters the 
cyclonic circulation flowing northward out of the South China Sea through the Formosa 
Strait. In winter, waters from the Philippine Sea enter the northern South China Sea through 
the Luzon Strait. A small amount of water flows northward, through the Formosa Strait, out 
of the South China Sea. Most parts of it flow westward and cyclonically from the Luzon 
Strait to Hainan and then it flows southward along the coast of Vietnam toward Singapore. 
It then flows eastward entering the Java Sea through the Gasper and Karimata Straits. There 
is another branch of water flowing southward along the coastlines of Luzon, Palawan, and 
Boreno to Natuna Island. On its way from Luzon to Natuna Island, waters continuously 
separate and flow westward spreading all over the central South China Sea toward the 


Asian coastline. 
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Figure 4A : Shallow Water Model Simulating the Summer Surface Circulation of the 


South China Sea (PohImann, 1987) 
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Figure 4B : Shallow Water Model Simulating the Winter Surface Circulation of the 


South China Sea (Pohlmann, 1987) 


This model predicts the seasonal varying surface currents quite reasonable when 
compared with the Wyrtki (1961) atlas. The currents are generally northward during winter 
and change to southward during summer. However, there is no jet branching during 
Summer, and the model simulated current pattern in the central South China Sea is different 
from the Naga Report observations (Figure 3A). Moreover, during winter the model 
simulated currents at the Luzon western coast, Formosa Strait, and the Dangerous Ground 
(Spratly Islands) are oppositely directed compared with the observations summarized in the 
Naga Report (Figure 3B). As pointed out by Pohlmann, the shallow water model might not 
be appropriate for the South China Sea due to the complex topography. 

A three-dimensional, primitive equation ocean circulation model with free surface 
developed at Princeton University (Blumberg and Mellor, 1987) will be used to study the 
South China Sea response to the monsoon system. The principle attributes of this model are 
as follows. 

*It contains an imbedded second moment turbulence closure sub-model (Mellor 
and Yamada, 1982) to provide vertical mixing coefficients. 

*]It 1s a Sigma coordinate model in that the vertical coordinate is scaled on the 
water column depth. 

*The horizontal grid uses curvilinear orthogonal coordinates and an “Arakawa C” 
differencing scheme. 

*The horizontal differencing 1s explicit whereas the vertical differencing is 
implicit. The latter eliminates time constraints for the vertical coordinate and permits the 
use of fine vertical resolution in the surface and bottom boundary layers. 

*The model has a free surface and a split time step. The external mode portion of 
the model is two-dimensional and uses a short time step (25 sec) based on the Courant- 
Fnedrichs-Levy (CFL) computational stability condition and the external wave speed. The 
internal mode is three-dimensional and uses a long time step (900 sec) based on the CFL 
condition and the internal wave speed. 


*Complete thermodynamics have been implemented. 


The grid size of this 23-level model is about 20 km in the horizontal. It will deal with 
both the shelf sea and the deep ocean regions in the modeling process and provide the 
numerical results to help us understand the dynamics of the South China Sea. The 
observations from the Naga Report are used for determining the open boundary conditions 
and for verifying the model results. The Levitus data set is used for initializing the model 
and for verifying the model results. 

The structure of this thesis is described as follows. Chapter II describes the 
background of the South China Sea which is based on the Naga Report and other 
supplements. Chapter IIT introduces the numerical model and the data sources used in this 
study. Chapter IV presents the model simulated results. Chapter V is the conclusions and 


the future recommendations. 
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Il. BACKGROUND 


A. CONFIGURATION OF THE SOUTH CHINA SEA 


The South China Sea is situated between the Pacific and Indian Oceans and formed by 
the Asian coastline and the chain of marginal islands of the Southeast Asian Waters. It 
extends about 3380 km from Taiwan to Singapore and about 1690 km from Hainan to 
Palawan. Table 1 shows the straits connecting the South China Sea with adjacent seas. It 
consists of many particular geographical features such as a deep sea basin, an ocean 
plateau, and a continental shelf, and these features make its regional bottom topography and 


oceanic appearances complex (Figure 2). 


TABLE 1: THE SRAITS CONNECTING THE SOUTH CHINA SEA WITH 
ADJACENT SEAS 


Adjacent Seas Si Depth) 


Se 
ete a Be 
ee ae | 


Its northern part is a continental shelf which includes the Gulf of Tonkin, Formosa 
















Strait, Hainan, and Taiwan. The average widths of this shelf are 152 km along mainland 


China and 10 km from most southern point of Taiwan. The Gulf of Tonkin deepens from 
the coast gradually to the center with a maximum depth of 70 m. On the other hand, the 
Formosa Strait is about 150 km wide, 60 m deep in the northern part and 50 m deep in the 
southern part (Gilg, 1970). 

The west and north central shelf areas are part of a submerged Asian continental 
margin which extend southward from the south of Formosa Strait and eastward from the 
coast of Vietnam into the central basin, respectively. It consists of a continental slope, 
numerous coral reefs and several ocean valleys. In the north central part there are numerous 
peaks of reef such as: Pratas Reefs, Helen Shoal, Paracel Island and Reefs, and 
Macclesfield Bank. The terrace east of Pratas lies between 2400 and 2950 m and is 
exceptionally deep. Most of the coral reefs occur on or within the terraces; however, 
Macclesfield Bank is an exceptionally large reef which is not related to a terrace. This bank 
and several large ridges to the south do not reach sea level and appear to extend with a steep 
slope to the basin floor on the east (Gilg, 1970). 

The central part is a deep sea central basin with a maximum depth more than 4300 m. 
The south central part is a plateau with reefs called the Spratly Islands or the Dangerous 
Ground. There are numerous reefs and islands on this plateau and they generally lie in 1470 
to 2010 m depth. Only the Reed Bank in the northeast is shallower and its base appears to 
lie in 300 to 500 m depth (Wyrtki, 1961). There is a channel with a sill depth about 2050 m 
through the reef region connecting the central basin and Palawan Trough. However, the 
Palawan Trough appears to connect with the central basin by way of a valley with a sill 
depth about 1400 m on the northeast side of Reed Bank, too. Moreover, an apron of 
sediments is present at the south end of the Dangerous Ground and coral reefs appear less 
numerous here (Gilg, 1970). 

The Gulf of Thailand is situated on the western side of the Sunda Shelf and its mouth 
opens about 338 km wide into the South China Sea. The distance from the head of the Gulf 
to its mouth is about 760 km with its greatest width about 500 km. Its bottom topography 


is shown in Figure 5A and its deepest place is about 80 m deep at the central part of Gulf. 


The northeastern coast is generally shallower and flatter than the southwestern coast. The 
deeper central gulf is separated from the South China Sea by two ridges. One of them 
extends southwestward at depths of less than 25 m for more than 100 km from Cape 
Camau. The other deeper ridge extends northeastward off Kota Bharu for a distance of 152 
km at depths less than 50 m. In the narrow deeper channel between the ridges a sill depth 
of 67 m was observed (Wyrtki, 1961). 

On the eastern side is a submerged island margin with the deep Palawan Trough of 
more than 3400 m depth. The slope of this margin is very steep near Palawan and Luzon. 
The average shelf widths of this margin are 10 km from Luzon, 50 km from Palawan, and 
84 km off Borneo (Gilg, 1970). It connects the Sulu Sea and the Pacific Ocean with the 
Balabac, Mindoro Straits and Luzon Strait, respectively. In the Luzon Strait, which is 
between Taiwan and Luzon, there is a depression in the form of a deep channel that extends 
northeastward across the submarine ridge called Bashi Channel (Figure 5B). This channel 
is 3100 to 3600 m deep throughout its length, except for sills of 1830 and 2380 m off the 
tip of Taiwan. Numerous small channels and valleys connect this deep channel with the 
Pacific Basin across the submarine ridge north of Luzon. One of them, northwest of Y’ Ami 
Island, has a sill depth of 2600 m and 1s probably the passage of main water interchange 
between the South China Sea and the Pacific Ocean. The other channels across this ridge 
have sill depths of 730 to 1300 m (Gilg, 1970). 

The southern part of the South China Sea is the Sunda Shelf which is situated between 
Malaysia, Sumatra, Java, and Borneo. It connects the Andaman Sea through the Strait of 
Malacca and connects the Java Sea through Karimata Strait and Gasper Strait, respectively. 
It is a very flat shallow trough with depths between 40 to 100 m at its central part and depths 
which decrease south of the equator to less than 40 m. However, there are several 
submerged river valleys which lie on the bottom of the Sunda Shelf between Borneo, 
Sumatra, and Singapore. They drain the land and join together flowing northward into the 


south central plateau area. These movements transport the riverine sediments and gradually 
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Figure 5A : The Bottom Topography of the Gulf of Thailand (Wyrtki, 1961) 
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Figure 5B : The Bottom Topography of the Luzon Strait (Shaw, 1989) 


make the water depth shallower in the southern part of the Sunda Shelf. Similar features 


are also found in the Gulf of Thailand (Wyrtki, 1961). 
B. ATMOSPHERIC FORCING 


1. Low-Level Atmospheric Circulation 

The South China Sea is an ideal region for monsoon development. Hess (1959) 
pointed out that the monsoon changes with season over Asia as a direct and simple response 
to the low continental temperatures in winter and the high continental temperatures in 
summer. Figure 6 shows the mean December and July 850-mb streamline analysis of the 
Southeast Asia monsoons (Cheang, 1987). There is a strong divergent anticyclone over 
Siberia (cold high) in winter and a convergent cyclone in summer. Therefore, the equatorial 
pressure trough crosses the equator twice each year corresponding to the north subtropical 
high and the monsoons over tropical low latitude regions will also change direction twice 
each year. 

In winter, northeasterly winds prevail over the whole South China Sea with an 
average magnitude of 9 m/s. In summer, southwesterly winds dominate most of the South 


China Sea and Gulf of Thailand with an average magnitude of 6 m/s (Wyrtki, 1961). 


2. Buoyancy Flux at Air-Sea Interface 
Rainfall plays a very important role in the South China Sea and the Gulf of 
Thailand. Figure 7 shows the average rainfall characteristics of the Southeast Asian Waters 
and Table 2 shows the yearly amounts of precipitation and the other heating properties, 
where Q,,,— Q, (annual total effective radiation energy - energy required for evaporation) 


is the available energy for heating processes. The units are: P in cm/yr, E in cm/yr, and 


Q.,— Q, in cal/cm*/day, respectively. 
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Figure 6 : The Mean 850-mb Streamline of the Southeast Asia Monsoons (a) 


December, (b) July (Cheang, 1987) 
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Figure 7: Average Rainfall (mm/yr) of the Southeast Asian Waters, Higher than 4000 


mm Dashed (Wyrtki. 1961) 
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TABLE 2: THE ANNUAL TOTAL PRECIPITATION AND HEATING 
PROPERTIES OF THE SOUTHEAST ASIAN WATERS (Wyrtki, 1961) 


Precipitation | Evaporation - 0.4-0, 


(cm/yr) (cm/yr) (cal/cm?/day) 


nll il el i 
Northern 
South China Sea 


Southern 100 
South China Sea 
~~ eee ee eee 





Heavy rainfall and its strong dependence on latitude are important features of the 
South China Sea. The precipitation rate decreases from 221 cm/yr in the southern South 
China Sea to 152 cm/yr in the northern South China Sea. On the other hand, the evaporation 
rate increases northward from 96 cm/yr in the southern South China Sea to 133 cm/yr in 
the northern South China sea. The northward increase in evaporation and decrease in 
precipitation leads to an heterogeneous fresh water influx into the South China Sea surface. 
The southern South China Sea receives more fresh water (106 cm/yr) than the northern 
South China Sea. Because of the heavy rainfall and large amounts of energy available for 


heating, the South China Sea and its adjacent seas are sources of warm, low salinity water 
(Wyrtki, 1961). According to Table 2, the annual total Q,,,—Q, (cal/cm?/day) values 


indicate that the energy available for heating is quite large in the southern South China Sea 


(100), the Strait of Malacca (127), and the Sulu Sea (107), and quite small in the northern 


South China Sea (38), Java Sea (77), and Andaman Sea (56). The north/south difference of 


the Q,,,— Q, values over the South China Sea is around 62 cal/cm?/day. 
C. OCEANIC CIRCULATIONS 


1. The Equatorial Current System 
In response to the monsoons and atmospheric features, the wind driven Pacific 
Ocean equatorial current system has at least four major currents. Three of them extend to 


the surface, the other remains subsurface. The three major upper layer currents are the 
westward-flowing North Equatorial Current between about 20 to 8 N, the westward South 


Equatorial Current from about 3N to 10S and the narrow North Equatorial 


Countercurrent flowing to the east between them. The fourth current, the Equatorial 
Undercurrent, flows to the east below the surface straddling the equator from about 2 N to 


2°S (Pickard and Emery, 1990). 


The North Equatorial Current (NEC) flows westward and is divided into two 


branches off Luzon at 13 to 14 N and they flow southward and northward, respectively 
(Pickard and Emery, 1990). The southward branch is the Mindanao Current. The northward 
branch belongs to the Kuroshio Current and flows along the coast of Luzon Island with a 


maximum velocity of 0.8 - 1.0 m/s, then enters Luzon Strait. After that it turns to the 
northwest and its main axis reaches as far west as 121 E with a maximum velocity of about 


1.5 m/s. There is one branch that often goes westward at 20 N and enters the South China 
Sea. However, most parts of it will return in the main axis of the Kuroshio Current (Nitani, 
1970). 

The NEC passes through the passages of the eastern submerged island margin of 
the South China Sea and mixes with the South China Sea water masses. Bashi Channel in 
the Luzon Strait is the only passage where the South China Sea waters can exchange with 
the lower layer and deep ocean water masses of the Pacific Ocean (Wyrtki, 1961). In 


general, the lower layer water flows in the opposite direction of the surface circulation. 
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Figure 8 shows the deep water connections between the Pacific and Indian Oceans and the 
flow of deep waters into the various basins. Deep water enters the South China Sea central 
deep basin through the Luzon Strait where it branches with one flowing northward into the 
Palawan Trough and the other flowing out of the South China Sea through the Mindoro 
Strait. There is no deep water connection between the Sulu Sea and the other deep basin 
waters. Therefore, the deep water masses in the South China Sea totally belong to the 


Pacific Ocean water system (Broecker et al., 1986). 


2. The Kuroshio Intrusion Through the Luzon Strait 

The seasonal variation of the Kuroshio intrusion through the Luzon Strait 
dominates the structure of water masses of the northern part of the South China Sea. Based 
on CTD data, Shaw (1989) hypothesized that the intrusion of water from the Kuroshio 
Current starts in late summer and reaches a maximum in winter. By late spring, the 
intrusion from the Kuroshio ceases, and the waters from the South China Sea begin to enter 
the region. Remnants of the intrusion water from the Kuroshio are trapped at the mouth of 
the Peng-Hu Channel. The water from the South China Sea eventually occupies the entire 
slope region south of the Peng-Hu Channel until the next intrusion of the Kuroshio begins 
in late summer. 

The intrusion of waters from the Kuroshio could reach at least the salinity 
minimum layer (about 500 m) and could extend as far north as the Peng-Hu Channel (Shaw, 
1989). Because the intrusion water is deep, the steep continental slope poses a strong 
dynamic constraint to the intrusion current. Consequently, the intrusion water can not 
continue through the Formosa Strait into the East China Sea. Recent hydrographic data 
showed that the shelf water in the Formosa Strait was separated from the deep water by a 
front south of the Peng-Hu Islands in winter. Consequently, it is unlikely that there is a 


branch of the Kuroshio intrusion in the Formosa Strait. 
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Figure 8 : The Deep Water Connections Between the Pacific and Indian Oceans 


(Broecker et al., 1986) 
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Many studies show that there are anticyclonic and cyclonic circulations on the 
northern South China Sea near the Luzon Strait (Figure 9). An analysis of the geostrophic 
current (Zhong, 1990) shows that the anticyclonic circulation at the continental slope 
consists of the South China Sea Warm Current and the intrusion water of the Kuroshio. A 
cyclonic eddy is induced periodically by the Kuroshio Current across the Luzon Strait (Liu, 
1992). Figure 9 also shows the seasonal variation of the surface current system of the Luzon 


Strait. 


3. The Formosa Strait 

The seasonal variation of the current in Formosa Strait is very evident. According 
to Fan (1982), the current is primarily composed of three different currents: Kuroshio 
intrusion, the China Coastal Current, and the South China Sea Current. The China Coastal 
Current originates near the coast of northern China and flows southward along the China 
coast. When the winter monsoon prevails, it is strengthened by the wind and flows all the 
way through Formosa Strait and encounters the Kuroshio intrusion at the Peng-Hu Islands. 

When the summer monsoon prevails, the South China Sea Current is strengthened 
and the Kuroshio intrusion is weakened by the southwest winds. Most of the South China 
Sea current flows through the Luzon Strait while a small amount enters the Formosa Strait 
flowing northwards. Shaw (1992) found that the cold shelf water in the East China Sea 
enters the Formosa Strait in fall and winter and also that the coastal discharge from the 


Yangtze River and rivers south of it produce a southward current into the Formosa Strait. 


4. The Java Sea 


The Java Sea is smaller than the South China Sea and the currents occupy its full 
width without the development of larger eddies and counter currents. From May to 


September the water flows to the west and from November to March to the east. In April 
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Figure 9 : Seasonal Variation of the Current System Near the Luzon Strait 


(Huang, 1988) 
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and October the direction of the flow reverses and eddies occur. Normally in these months 
a current towards the east prevails off the coast of Borneo. In the small passages between 
Borneo and Sumatra, the winds are fresh and the speeds of the currents often reach 1 m/s, 


especially during the full north monsoon (Wyrtki,1961). 


5. The Strait of Malacca 
Through the Strait of Malacca a small exchange of water is possible between the 
Sunda Shelf and the Indian Ocean. The water movements are in general directed towards 
the Indian Ocean and are strongly related to the surface gradient of the sea level through 
this strait. In the Strait of Malacca the period of stronger flow is from January to April, 
during the northeast monsoon, but it is chiefly caused by the low sea level in the Andaman 


Sea in this season (Wyrtki, 1961). 


6. The Gulf of Thailand 
Very few observations have been made in the Gulf of Thailand but from available 
data (Wyrtki, 1961) water movements seem to be weak. During the summer monsoon 
season they are anticyclonic and when the winter monsoon prevails they are cyclonic. 
Neelsasri (1979) found that density forcing is as important in driving the circulation as the 
mechanical forcing of a uniform wind field. He also reports that the currents of the South 
China Sea induce a shear flow near the open side of the Gulf and that they are confined to 


the outer half of the Gulf. 


7. The South China Sea 
The circulation of the South China Sea is mainly determined by the monsoons and 
the configuration of the land masses surrounding the sea. The ship drift currents of the 
summer and winter monsoon seasons are shown in Figures 3A and 3B (Wyrtki, 1961), 
respectively, and summarized as follows. 
In June, the largest part of the South China Sea water masses pass south of Taiwan 


into the root of the Kuroshio; a smaller part flows northwards through the Formosa Strait 
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joining the Kuroshio Current. Because the supply of water to this current from the Java Sea 
seems to be insufficient (Wyrtki, 1961), a current develops off the coast of Borneo towards 
the southwest with a current speed between 0.12 m/s and 0.25 m/s, which turns back to the 
north in the region of Natuna Island. The water of this current comes partly out of the Sulu 
Sea, and partly from a counter current, which is formed in the central parts of the South 
China Sea. This current must develop in close connection with the westward intensification 
of the main current and be caused dynamically (Wyrtki, 1961). Along the western coast of 
Luzon water flows northward from the Sulu Sea with a surface current speed of between 
0.12 m/s and 0.25 m/s. The surface current speed is more than 1.0 m/s off the coast of 
Vietnam (Wyrtki, 1961) and about 0.5 m/s off the Hainan coast. In the Gasper and Karimata 
Straits, the surface current speed is about 0.75 m/s and the current speed of the jet current 
which flows from Vietnam to the Luzon Strait is about 0.25 m/s to 0.30 m/s during this 
season. 

In December, the northeast monsoon 1s fully developed over the South China Sea 
and the currents are strong. Waters from the Kuroshio Current and the East China Sea enter 
the South China Sea flowing southwestward along the coastline of Asia and flow out of the 
South China Sea through the Gasper and Karimata Straits. The surface current speeds are 
about 0.50 m/s off the coast of Hainan, about 0.75 m/s off the coast of Vietnam, and about 
0.30 m/s in the Gasper Strait. The westward intensification is stronger than during the 
southwest monsoon, when a deflection of the current occurs to the east caused by the wind. 
The counter current in the center of the South China Sea 1s present and lies farther to the 
east because of the strength of the main current. From the Sulu Sea water flows into the 
South China Sea and is mixed into the main current. A current boundary exists between 
Vietnam and Borneo and a very weak (less than 0.12 m/s) southwestward flow occurs along 
the coast of Borneo which encounter the waters from the north flowing northward along the 


right hand side of the current boundary. 


26 


D. WATER MASSEES 


1. Mass Transport 

Wyrtki (1961) showed that on the western side of the Luzon Strait the water 
masses can be divided into the warm, less saline surface layer water, the high salinity (34.8 
-35.2 psu) Subtropical Lower Water, the low salinity (34.1-34.5 psu) Intermediate Waters, 
and the original Pacific Deep Water (34.65-34.7 psu). The salinity maxima of the 
Subtropical Lower Water is about 150 m deep and the salinity minima of the Intermediate 
Waters is about 400 to 500 m deep. According to Wyrtki (1961), in the Luzon Strait these 
two water masses move in opposite directions and the depth of this inversion is about 200 
m. The center of this inversion lies relatively uniformly at a temperature 16°C. 

Wyrtki (1961) also found that during the winter monsoon season surface water 
from the Pacific Ocean is transported into the South China Sea through the Luzon Strait. 
The vertical extent of this layer is about 200 m. This Subtropical Lower Water spreads 
chiefly along the western side of the South China Sea following the surface currents. Parts 
of the Intermediate Water flow out of the South China Sea in depths between 400 and 900 
m during this season. During the summer monsoon season the conditions are reversed. The 
surface water flows out of the South China Sea and below about 300 to 400 m water masses 
flow into the South China Sea. 

According to Nitani (1970) the North Pacific Intermediate Water forms the lower 


part of the Kuroshio Current. The transport of the Kuroshio Current off Luzon is about two- 


thirds to three-fourths that of the Mindanao Current, which is about 30-45 Sv (10°m?/s), 
and is smaller than that of the Kuroshio Current off the south coast of Japan. Figure 10 
shows the mass transport in June and December of the Southeast Asian Waters and Table 
3 shows values of mass transports for different currents. According to Table 3, the mass 
transport of the Kuroshio Current to the east of the Luzon Strait is the largest transport in 
the area. However, in the South China Sea the largest transport occurs off the Vietnam 


coast. 


Dad 





Figure 10A : Summer (June) Mass Transport (Sv) of the Southeast Asian Waters 


(Wyrtki, 1961) 
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Figure 10B : Winter (December) Mass Transport (Sv) of the Southeast Asian Waters 


(Wyrtki, 1961) 
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TABLE 3: MASS TRANSPORT (Sv) OF CURRENTS IN THE SOUTHEAST 
ASIAN WATERS (Wyrtki, 1961) 


[Currents | Feb. | Apr | gun] Aug | Oct | Dec 

Kuroshio Current 24.5 SHS, 35.0 30.5 2740 215 
(Northwards positive) 

Off Vietnam Current 5.0 1) -3.5 -3.0 2.0 5.0 
(Southwards Positive) 

Timor Current 1.0 15 1.0 1.5 les L.5 
(Westwards positive) 

Luzon Strait Current 2:5 -3.0 -2.5 0.5 3.0 
(Westwards Positive) 












2. Temperature and Salinity 

According to Liahude (1979) during the winter monsoon season the northeast 
monsoon brings cool (about 26°C), high salinity (about 33.5 psu) water from the northern 
South China Sea into the southern South China Sea. During the summer monsoon season 
the southwest monsoon brings warm (about 29°C), low salinity (32.5 psu) water from the 
Java Sea into the southern South China Sea. The waters in the southern South China Sea 
and the Java Sea are warm and less saline. Thus, the variations of water properties in the 
South China Sea are mainly determined by the cool, high salinity waters from the north in 
the winter monsoon season. On the other hand, the river runoff plays an important role in 
determining the surface salinity in the Gulf of Thailand, the Gulf of Tonkin, and the regions 
around Borneo (Wyrtki, 1961). Large salinity gradients also exist in these regions. 

The Levitus (1982) surface temperature distributions in the South China Sea for 
summer and winter are shown in Figures 11A (June), 11B (December), and seasonal 


surface salinity in Figures 12A (June), and 12B (December). 
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Figure 11A : Levitus Surface Temperature Distribution (June); Interval = 1°C 
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INTERPOLATED LEVITUS TEMPERATURE CONTOUR 
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Figure 11B : Levitus Surface Temperature Distribution (December); Interval = 1°C 


32 


AGNLILV 1 


INTERPOLATED LEVITUS SALINITY CONTOUR 












MIN. MAX. INTERV. (ppt) MONTH 6 
31.50 34.50 0.25 Depth = Om 
100E 105E 110E 115E 120£ 

25N r— . : . 25N 








100E 105E NOE 115E 
LONGITUDE 


Figure 12A : Levitus Surface Salinity Distribution (June); Interval = 0.25 psu 
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INTERPOLATED LEVITUS SALINITY CONTOUR 
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Figure 12B : Levitus Surface Salinity Distribution (December); Interval = 0.25 psu 
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Summer: In June (Figure 11A) surface temperature is homogeneous (about 
29°C) in most part of the South China Sea, except in the Gulf of Tonkin and the Formosa 
Strait. In the Formosa Strait they are between 26.5°C and 28°C. 

The surface salinity 1s between 33.0 psu and 34.0 psu in most parts of the South 
China Sea (Figure 12A). Exceptions occur in the Formosa Strait, the Gulf of Tonkin, the 
Gulf of Thailand, and the coast of Borneo. Surface salinities are about 30.5 psu to 32.5 psu 
in the Gulf of Thailand and 31.5 psu to 33.25 psu in the Gulf of Tonkin, respectively. In the 
Formosa Strait surface salinity is more than 34.0 psu. Off the coast of Borneo surface 
salinity is between 32.5 psu and 32.75 psu. 

Winter: In December (Figure 11B) the isotherms spread southwestward 
increasing from 22°C in the Formosa Strait to 28°C off the coast of Borneo. The largest 
surface temperature gradients occur in the Formosa Strait (20.5°C to 22.5°C). In the Gulf 
of Thailand the surface temperature is higher than 27.5°C. In the Gulf of Tonkin, however, 
surface temperatures are between 23.0°C and 25.5°C. 

The surface salinity increases with latitude in winter from 33.0 psu in the southern 
South China Sea to 34.0 psu in its northern reaches (Figure 12B). Large salinity gradients 
occur in the Gulf of Tonkin, the Gulf of Thailand and the western side of Borneo. In the 
northern South China Sea, surface salinity 1s between 33.5 psu and 34.0 psu. In the Gulf of 
Tonkin and the Gulf of Thailand, surface salinity varies from about 31.0 psu to 33.5 psu. In 
the southern and central South China Sea, surface salinities are homogeneous between 
33.25 psu and 33.5 psu. The surface salinity is more than 33.5 psu in the Gasper and 


Karimata Straits. 
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Il. NUMERICAL MODEL 


A. INTRODUCTION 

The Princeton turbulent closure model (Blumberg and Mellor, 1987) was 
implemented for the South China Sea at the University of Miami and the Naval 
Postgraduate School. The model grid is shown in Figure 11. The model domain is 99 E - 
121 E, 3 S - 25 N which includes the South China Sea and the Gulf of Thailand. With its 
20 km horizontal resolution and 23 vertical sigma coordinate levels, it has 125 x 162 x 23 
grid points. By using the mode splitting technique and finite difference formulation 
introduced by Blumberg and Mellor (1987), the governing equations for momentum, heat, 
and salt can be simplified to predict four different prognostic variables: surface elevation, 
1, temperature, 7, salinity, S, and velocity, (U,V,W). 

The data sets needed in this simulation are bathymetry, surface wind stress, 
temperature, salinity, and mass transport at the lateral open boundaries. NAVOCEANO 
dbdb5 bathymetry, Hellerman wind stress, and Levitus temperature/salinity data sets are 


available for this modeling research. The original data resolutions are as follows. 

* NAVOCEANO dbdbS5 Bathymetry (5 x5) 

* Monthly Levitus temperature/salinity (1 x 1 ; 19 levels). 

* Annual Levitus temperature/salinity (1 x 1 ; 33 levels). 

* Monthly Hellerman wind stress (2 x 2 ). 

All the data used in this simulation were interpolated into the model grids. 

Since the Levitus monthly data sets only have 19-levels, the data below the 19th-level 
(levels 20-33) were taken from the mean annual climatological Levitus data. Here, a data 
smoothing technique was used to prevent an unrealistic jump between the 19th and the 


lower levels. 


36 


LATITUDE (N) 


2 eeeeno ae 









Soa OSes sae 
« S22 e2eeeees ca 
222 e2e eee Labed 
SSeauo a! 
a a 

SS eS 
2oe@ Sees Seasee 
sf2enee0@ 










rid 
1 SaaS e228 868808 
OT hed dedende-chrah lyadedendndhadhadpadhoad 
= 2S SSeS eaeen 
+7 SR OSSSPSTSSCES 
oP MP ORSGSSSSSOSOOCAS TTS SSSSSHRSE ASSESSEPSSSSS@ASSOCSS 
o 


= * 
2ee8c8. Gee) Yt ewes Oeeeee 

dat ee 0Ult(‘<i‘éwmR err ttt PT LTT Tt Te tt 

a@e Sega 2228 O28 SOSEGEo Feo BSS OD 
€ beer Oe OSGeo 8B SSS SSS S28 2 PSF 2 SSS 2S SS SS SSCS SSS SE SSFSS5 SVOSCSHOSISVOSSATS SOS. 
tritttotiietttm oiuettrmsesrittitittiititrirrcritcitctttittttttt Tritt tT Tet er 
SP SSGeeehsGnoseoon..€§ Be" * ®eeenrea e 
OPS OSOSREOSSSSSGE. ” TEaaeSorheae Seacou 
Seene et eSeeee ee ene. 6 6lvt ewe aereen - sea 
29 8S SSeS SSE baer SSSSoPeaaeg aoe cane 
[SL ddetetietinttentttintinteateee pte Pi 

S2oone 


S29 seers eccte ns! ii ee jo ocS ceases ee oem 
Se SSCS OSCSTsS= SS eeeeaSeeoqg nae 
seaee@ene@eweasae  - -ceecece Geecenee 
s#@ESS 288806890280 






sacceaeoeasses BE “eee sec ees eaeauvess 
eseeeuesececorus! EE *¢ Seu even e209 eeER SER 
© OS 8S 6 8608 OO eS 8 28 6 O86 O66 8 O00 FOES OE BS OF OF OS SSE OO 64688 SSE FF SS SSE ES S8286; 
Tiitittdidde ttre  fereiietiss ts. ti tte 2 tt 
SES SSSA SsOwues5 FePG Fevers SOF SSOSSSOSESesleaoe 
PT Lt dpb dlp dehedebedadededudLed-debedehepet doch. debeddebotebedohebebededed 
SF CSAS SPS SPR OOO RSD SSSPS SAS OSS SSS SSL SG Ss SSP Ss SPF Ve st@c SSP Ss OES SIO BS SCTS# SSeS SOee sow 
PSs SS SSSE SSSR SRE SOF CSSSESSHOSSESSSESTIES 
SPSS SC OS CF OS OS OSS 6 OF SFOS SES 28SEC OBOE SBS BE ESS 46 OOS O26 2ESCOSHSSSESOSVSOSEEBDSOE 
_ DUST PS PAFETISSSSASSHRISOA SSS SRI SOSAASSIVSSVSISCSeue SQeracossssoa 
seeanecoaenese! 





—- 
On 


See 8 SRF 6 SSS tS ee Se 
aS OPP Aaa SseHSaeses 
SSeS S test eeres eae 
oa 2 SF O® OWS Gees & & a= 
PO BSSSS2 CS SSESESS SSSPa 
6 F228 2 SSS PSST] STP OMS SS Baw. 
029 892280 2 CB 2SS2SE 29 S248 280 


=i 
! 
; 
| 


Ga Ga eae ep Gh eis mh eres Ss © os ea GO ee ge” 
~ eS ES See SS eee eee ee | 

1889 3SS 65592 5882 PSS SSES6ERR 
see * SOs aso mag aa a 
ASS FPeStesoS S2SG20 Soom a 
PISMO SveSCsgSsesosesssscga 
RP SAA SSPOs SMSwre Se Sees oemraaGa and @ 2 op ab GD a a> Sate ee ee 

SOOSE SSE SS COS SSS FESTESS G2 SSTSS 82880 


SSS Ee Se GSS + FOS OS & ewe © © GS © ee eS © a ee 
SSS SSS S SH SE SSSESIESSECSES FOSTS“ES OCB 
<= © SS SES SFLS SCSSSelOS SSeS SeOessesaeeeo 

_ 8 PS SES SS STS SSSSS SESS SES 
SO © @ SSSSSSSESSOSOSTSOEESSEOD 


105 —i(<ié‘ ONC 115 420. 
LONGITUDE (E) 


Figure 13 : Model Grid 
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Hellerman and Rosenstein (1983) used a quadratic bulk aerodynamic formulation to 
calculate the surface wind stresses. Over 35 million oceanic surface observations exist, 


extending from 1870 to 1976 and archived by the National Climatic Center in Asheville, 


NC, designated TDF-11, and grouped into data set with 2 x 2. resolution. 


B. MODEL FEATURES 


The model features are listed as the follows: 

* Primitive equation model. 

* Second order turbulent closure. 

* o coordinate (vertical). 

* Curvilinear orthogonal coordinates (horizontal). 

* Free Suthaces 

* Split time steps (internal mode = 900 sec.; external mode = 25 sec.) 
* No slip hard boundary conditions. 


* No fluxes on rigid boundaries (bottom, side walls). 


C. MODEL GOVERNING EQUATIONS 


According to Blumberg and Mellor (1987), this model uses the implicit numerical 
scheme in the vertical direction and a mode splitting technique in time which has been 
adopted for computational efficiency and which separates the external mode (vertically 
integrated equations) from the internal mode (vertical structure equations). Its numerics 
will accommodate the highly time-dependent and often nonlinear processes of coastal 
upwelling and eddy dynamics. By using the hydro-statistic and Boussinesq assumptions 
and considering a system of orthogonal Cartesian coordinates with x increasing eastward, 
y increasing northward, and z increasing upward, the governing equations of this model, 
which form the basis of the circulation model, describe the velocity and surface elevation 


fields, and the salinity and temperature fields. 
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The free surface is located at z = N(x, y,f) and the bottom boundary is at 


z = —H (x,y). The V is the horizontal velocity vector with components (U, V) and V is 


the horizontal gradient operator. 


1. The Continuity Equation 


Eye 0 (1) 
Oz 
2. The Momentum Equations 

dU = dU ete laa dU 
Ry Me ed —-fV = Ooener | ac (2a) 

ov > OV oe clare OV 

pgs = — 

Oz (2c) 


with p, the reference density, p the in situ density, g the gravitational acceleration, P the 
pressure, K,, the vertical eddy diffusivity of turbulent momentum mixing, and f the 
Coriolis parameter.The pressure at depth z is: 

PG yc, = FP + SP te 10 Oy, 2 dc 


4 ° 
where 77 1s the free surface and the atmosphere pressure, P.,,,,, 1S assumed constant. 


atm? 


3. The Conservation Equations for Temperature and Salinity 


OF net We = 2K) +F 3 
Ee Dz oz °8¥ag) the fale 
ase aS) Pees 

yu a Wa = 5, (Kus-) + Fs (3b) 


where 9 is the potential temperature (or in situ temperature for shallow water applications) 
and S is the salinity. The density is computed by using the equation of state of the form 
0 = p,(8, 5) where p, is the potential density. K,, denotes the horizontal eddy diffusivity 


for turbulent mixing of heat and salt. 
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All of the motions induced by small-scale processes not directly resolved by the model 
grid are parameterized in terms of horizontal mixing processes. These unsolved processes 


are denoted by F terms in the following equations 


_ 3 QU |. oV 
Lee = oe Cy aha ai ve 7 oe ayia Be) (4a) 
A eV. Q oU | aoV 
oy = Seat Ra a) (4b) 
i 0 0 0 
Fos = a, Ana, (8; S)) Tay lari S)) + Res (4c) 


where A, and A, are the horizontal diffusivities, which are the horizontal turbulent 
exchange coefficients for momentum and heat/salt mixing and the atmospheric 
thermodynamic forcing, respectively. Rg . is the term representing relaxation to the Levitus 
temperature and salinity. Although the F terms allow variable coefficients, their values as 
low as 10(m’/s) have been used successfully in various applications. We use the 


Smagorinsky diffusivity for ype eee and its formula fe, 1 


uv = 0.1 (AxAy) a4 ) 405 (+0) +) ) 


D. TURBULENCE CLOSURE 
The vertical mixing coefficients, K, and K, in (5a),(5b),(6a),(7a) are determined by 


the second order turbulence closure scheme (Blumberg and Mellor, 1987). According to 


(Sa) and (5b), they will be characterized by equations for turbulent kinetic energy (TKE), 
q’/2, and the turbulence macroscale (turbulence mixing length), /. 


aq? og -2 7 24? 3U.2_aVv.? 





or 
28 Kx - +r, (Ja) 
aig +V- Vig? + wella = 2 (KS) 1B K, (LY) +X) 9+ 
ae ks - Z W+F, (Sb) 


se 2 
wmere W = 1+ 6, () is the wall proximity function and L_ can be calculated by 


L" = (n-z)7!+ (H+z)7!. Near the surfaces the Von Karman constant k = 0.4, and 


thus //K=L,W = 1+ E,, and far from the surface / « L, W =1.The first term on the right 
hand side of (5a),(Sb) 1s diffusion, the next two terms are shear production and horizontal 
mixing. The horizontal mixing terms, F, and F,, can be parameterized analogously to 
temperature and salinity by using (4a),(4b). The vertical mixing coefficients (K,,, K,,, and 
K,) are determined from Ky = /qS,,, Ky = 1qS,, K, = lqS,, where S,,, S,, and S, are 

dU a __ Op 


stability functions that are analytically derived and depend on variables ee and 55° 
&j 


4; Po,» and lf. 


E. SURFACE AND BOTTOM BOUNDARY CONDITIONS 


The boundary conditions at the free surface (z = 1 (x, y) ) are: 


Ku (ger ge) = (Top ty) (6a) 
a 
q = BiUis (6b) 
Gi =O (6c) 
w= UStevet a (6d) 
where (T,,, T,,) 18 the surface wind stress vector and U,, is the magnitude of the ocean 


surface friction velocity. B is an empirical constant arising from the turbulence closure 
relations. On the bottom and the side walls of the basin, the normal gradients of 8 and S$ 
are zero. No energy flux could pass through these boundaries. The boundary conditions at 


the bottom (z = —H (x, y) ) are: 


4] 


dU aV 


PoKu Ee ) aD = ge t,,) (7a) 
2 
g = BiUi, (7b) 
gi=0 (7c) 
oH oH 
W, = (-U,) 5. - Vis (7d) 


where the —H (x, y) is the bottom topography and U,, is the bottom friction velocity 
associate with the bottom frictional stresses (T, ,T fe According to the logarithmic law, the 
bottom stresses can be calculated by the equation: 
ig = p,CrIV, V, 
The drag coefficient, Cp, is determined by the equation: 


| In(H+z,).~ 


( Gage 
Sik iaaeay 
where z, = 0.01m and Cp = 0.0025 are recommended by Blumberg and Mellor (1987) 


and will be used in this simulation. z, and VA are the grid point and the corresponding 
velocity in the grid point nearest the bottom, respectively. If enough resolution data is 
provided and associated with the turbulent-closure-derived K,,, the velocity in the lower 


boundary region can be calculated by the equation: 


- 
V= Z Jo ES : 
(a i 


The value of roughness length, z,, depends on the local bottom roughness. 





F. o COORDINATE SYSTEM 
In this simulation the sigma coordinate system is used. The governing external and 
internal mode equations are transformed from (x,y,z,t) to (x, y, G, f)coordinates, 


z—T 


ei =f 
H+Mn 





xX =xX,y¥ = y, G— 
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and the o ranges from o =O at z = 1 (surface) to o =-1 at z = —H (bottom). Thus, all 
the derivatives in the dynamic equations will be transformed from the Cartesian coordinate 
system to the new system and used in this modeling research. By using a cubic spline 
interpolation, the sigma coordinate values and the Cartesian coordinate values can be 


transformed to each other as needed. 


G. INITIAL CONDITIONS 
The model is integrated from 30 April (t=O). The initial conditions of the model are 
U=0, V=0,W=0 
T = Theo 
S = Sho 
where 7],-9, Sl, are the temperature and salinity on 30 April, which are obtained by 


interpolating the Levitus monthly data. 


H. OPEN BOUNDARY CONDITIONS 

A lateral open boundary is any place where the water can flow through. The mass 
transports through the Strait of Malacca, Banka Strait, Balabac Strait, and Mindoro Strait 
are quite weak. Therefore, for simplicity, these straits are kept closed in the model 
simulation. The three other boundaries are remain open, the Luzon Strait, Formosa Strait, 
Karimata and Gasper Straits. The water mass properties (temperature and salinity) are 
given according to the monthly Levitus data. The bi-monthly specified mass transports are 


given in Table 4. 


TABLE 4: MONTHLY MASS TRANSPORT (Sv) AT THE OPEN BOUNDARIES 
Gasper&Karimata Straits 4.4 -4.0 -3.0 1.0 4,3 
(Eastwards Positive) 
Luzon Strait -3.5 3.0 Ze -3.4 
(Eastwards Positive) 
Formosa Strait Lo 0.5 -0.4 
(Northwards Positive) 
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I. ATMOSPHERIC FORCING 

The atmospheric forcing includes wind forcing and thermodynamic forcing. In this 
simulation the Hellerman wind stress has been used as the wind forcing data. Figure 14 
shows the Hellerman wind stress distributions over the South China Sea. In summer, winds 
from the Indian Ocean and Australia meet at Natuna Island blowing northward and 
counterclockwise all over the South China Sea. Winds originating from Australia are very 
weak and extend only to Natuna Island. The air stream from the Indian Ocean 1s stronger 
and dominates the atmospheric circulation during this season. The strength of winds in May 
(Figure 14A) is about five times weaker than in June (Figure 14B). In May, northerly winds 
can reach up to 20°N where they converge with the southwestward winds from the East 
China Sea and blow westward from Luzon to Hainan. In June, the prevailing winds over 
the South China Sea blow northeastward and counterclockwise from Vietnam to the Luzon 
Strait. In winter, winds from the Pacific Ocean blow southwestward from Taiwan to 
Singapore and then turn left toward the Java Sea. In November (Figure 14C), winds are as 
strong as in December (14D), but their southerly extent is limited to 10°N whereas in 
December they extend to Natuna Island (5°N). 

Heat fluxes are very important for the equatorial region. According to Levitus (1982), 
the South China Sea has an annual mean temperature of 28°C, higher than its adjacent 
Pacific Ocean waters. Therefore, the South China Sea water is an extensive heat source for 
the adjacent regions. In this simulation, the relaxation to the Levitus temperature and 
salinity data were included in the thermodynamic forcing. This prevents any extreme 
deviation from climatological mean values and guarantees that the South China Sea will act 


as a heat source. 
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Figure 14D : Winter (December) Hellerman Wind Stress Distribution 
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IV. MODEL RESULTS 


In this study the Princeton South China Sea (PSCS) model is used to simulate the 
seasonal circulation in the South China Sea. The Naga Report (NR) observations (Wyrtki, 
1961) and the Pohlmann shallow-water (PSW) model simulated results (Pohlmann, 1987) 
are used to compare the performance of the PSCS model in the South China Sea. Because 
the PSW model provided the simulated summer and winter mean surface circulation 
patterns and the NR observations provided even-monthly surface circulation depictions, in 
this study the surface circulation patterns of these two sources are compared to verify the 
PSCS model results quantitatively. The PSCS model simulated surface temperature and 
salinity are also discussed to help understand the dynamics of the surface circulation in the 
South China Sea. Moreover, the comparisons of surface temperature and salinity between 
the PSCS model and Levitus data (Levitus, 1982) are examined to determine the deviation 
of the model results from the observations. 

The PSCS model started running on 30 April, a time the southwest monsoon is very 
weak, and was run for one year. The spin-up time of the PSCS model is about 60 days, 
longer than the 15 days required for the PSW model. The longer spin-up time is due to the 
fact that the PSCS model’s horizontal and vertical resolutions (20 km, 23 levels) are two 


times higher than the PSW model (50 km, 12 levels). 


49 


A. SURFACE CIRCULATION 


1. PSCS Model Results 

Summer: Because the winter monsoon changes to the summer monsoon in the 
early part of summer, the time of 10 May may be considered a transition period with 
monsoon conditions being quite weak (Figure 14A). In response to the weak wind stress on 
10 May, the surface current speed is generally quite low over the whole area (Figure 15A). 
The maximum current speed, about 0.2 m/s, is found in the Gasper and Karimata Straits 
where waters from the Java Sea enter the South China Sea. On 30 June, summer monsoon 
prevails (Figure 14B) and surface current speed becomes much faster in the whole region 
(Figure 15B). Waters from the Java Sea enter the South China Sea through these two straits 
with a current speed of about 0.7 m/s. Figures 15A and 15B both show that the surface 
circulation in the Gulf of Thailand is anticyclonic in the summer monsoon season. On 30 
June, the waters from the Java Sea divide into two currents between Borneo and Singapore. 


The largest amount of water flows northward and separates into three branches at about 


109 E, 16N flowing northward along the Vietnam coast, northeastward toward Luzon, 
and southwestward toward Natuna Island. The northward branch flows into the Formosa 
Strait. The second branch appears as a strong jet flowing towards Luzon entering the root 


of the Kuroshio Current through the Luzon Strait. A third branch meets the Borneo Coastal 


Current at 5N, 114E flowing northwards and combines with the jet current from the 
second branch. The current speed is more than 1 m/s off the coast of Vietnam, about 0.5 m/ 
s off the eastern coast of Hainan, and about 0.2 m/s for the jet current from Vietnam to the 
Luzon Strait. 

Winter: Because the summer monsoon changes to the winter monsoon in late 


fall, the winter monsoon is weaker on 5 November (Figure 14C) in comparison to the 
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Figure 15A :Simulated Summer (10 May) Surface Circulation 
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Figure 15B :Simulated Summer (30 June) Surface Circulation 
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winter monsoon on 10 December (Figure 14D). In response to monsoons on 5 November 
and 10 December, the surface circulation patterns are similar, but the surface current speed 
is lower on 5 November (Figure 15C) than on 10 December (Figure 15D). Figures 15C and 
15D both show that surface circulation in the Gulf of Thailand is cyclonic in the winter 
monsoon season. On 10 December, waters from the East China Sea enter the South China 
Sea through the Formosa Strait. Small branches of water from the Kuroshio Current enter 
the South China Sea through the Luzon Strait. The water masses from the East China Sea 


and Kuroshio Current flow southwards along the Chinese and Vietnamese coastlines. This 
current passes Singapore at about 2 N, 104 E and flows southeastward entering the Java 


Sea. The waters around Palawan separate into two branches at about 114 E, 10 N. One 


branch turns right flowing northward and cyclonically along the western coast of Luzon. 


Then, it turns left at 18 N, 116 E and flows westward towards the coast of China where it 
joins the main stream of the coastal current flowing southwards along the Hainan and 
Vietnam coastlines. The second branch flows westward spreading cyclonically around 
Natuna Island, passing Singapore, and flowing out of the South China Sea through the 
Karimata and Gasper Straits. The surface current speed is about 0.5 m/s to 0.8 m/s off the 


eastern coast of Hainan and 0.7 m/s off the coast of Vietnam. 


2. Comparisons 
The surface circulation pattern comparison among the NR observations, the PS W 
model, and the PSCS model is investigated to verify the PSCS model and to examine the 
relative performances of the PSCS and PSW models. Tables 5 and 6 show the comparisons 
of summer and winter surface circulation patterns among the PSCS model, the PSW model, 
and the NR observations. Although the PSCS model exhibits some disagreements with the 


NR observations in the central part of the South China Sea and the region around Palawan 
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TABLE 5: COMPARISON BETWEEN THE NAGA REPORT AND MODELING 
RESULTS (SUMMER) 


Naga Report 
(Wyrtki, 1961) 


Locations 
Formosa Strait 


Luzon Strait 


Western Coast of 
Luzon 


Coast of Vietnam 


(June) 


Out Flow: The south 
China Sea warm cur- 
rent flows into the 
Strait. 








Out Flow: Currents 
flow along the western 
coast of Luzon, cross 

the southern tip of Tai- 
wan and then enter the 
Kuroshio main stream. 


Currents flow north- 
wards along the west- 
erm coast of Luzon. 


(1).Longshore currents 
flow northeastward. 
(2). A branch of cur- 
rents turn backward 
from off Vietnam to the 
south. 

(3). A jet current flow 
toward the Luzon 
Strait. 


Region around Palawan | Waters from the Sulu 
Sea spread into the 
Dangerous Ground. 


Central South China Currents flow north- 
Sea eastward from the 


coast of Vietnam. 
















Shallow-Water 
Model(Pohlmann) 
(Summer) 


Out Flow: The South 
China Sea warm cur- 
rent and the branch of 
Kuroshio current enter 
the Strait. 


In Flow: A Kuroshio 
branch flows into the 
South China Sea 
through the northern 
Luzon Strait. 

Out Flow: The South 
China Sea warm cur- 
rent flows out of the 
South China Sea 
through the southern 
Luzon Strait. 


Currents flow south- 
wards along the west- 
ern coast of Luzon. 


The currents spread all 
over the central South 
China Sea and flow to 
the western coast of 
Luzon. 


Waters from the west- 
ern side flow towards 
Palawan. 


Currents flow east- 
wards from the coast of 
Vietnam. 









Princeton South 
China Sea Model 
(30 June) 


Same as the Naga 
Report. 












Same as the Naga 
Report 















Same as the Naga 
Report 


Same as the Naga 
Report. 


Waters from the west- 
ern side flow towards 
Palawan. 


Same as the Naga 
Report. 


Karimata Strait In Flow Same as the Naga Same as the Naga 
Report. Report 

Gulf of Tonkin clockwise circulation (not clear) Same as the Naga 
Report. 

Gulf of Thailand clockwise circulation Same as the Naga Same as the Naga 
Report Report. 
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TABLE 6: COMPARISION BETWEEN THE NAGA REPORT AND MODELING 
RESULTS (WINTER) 


Shallow Water 






















Princeton South 
Model (Pohlmann) China Sea Model 
(Winter) (10 December) 
Out Flow: A branch Same as the Naga 


of Kuroshio current Report. 
enters the Strait. 


Naga Report 
(Wyrtki, 1961) 
(December) 
In Flow: Waters from 
the East China Sea flow 


into the South China 
Sea. 


Locations 
















Formosa Strait 













Luzon Strait In Flow: A small In Flow: Large Same as the Naga 
branch of Kuroshio cur- | amount of Kuroshio Report 
rent enters the South current flows into the 


China Sea South China Sea. 


Western Coast of Current flow north- Currents flow south- 
Luzon wards along the west- wards along the west- 
em coast of Luzon. ern coast of Luzon. 


Coast of Vietnam Longshore currents Same as the Naga Same as the Naga 
flow southwards. Report. Report. 

Region around Palawan | Waters from the Sulu Waters from the north Same as the Naga 
Sea enter the South flow southward along Report 


China Sea flowing the Palawan coastline. 
northward. 

















Same as the Naga 
Report 



















Currents flow west- Waters flow north- 

wards from the western | wards, turn left to the 

coast of Luzon and coast of Hainan and 

spread all over the cen- | then flow to the coast of 

tral South China Sea. Vietnam. Then it enters 
the main stream go 
southwards. Unlike the 
Naga Report the waters 
continue southward 
and do not form a cur- 
rent boundary. 


Karimata Strait In Flow Same as the Naga Same as the Naga 
Report. Report 

Gulf of Tonkin counterclockwise circu- | (not clear) Same as the Naga 
lation Report. 

Gulf of Thailand counterclockwise circu- | Waters inflow from the | Same as the Naga 
lation South China Sea Report. 


Currents flow counter- 
clockwise and current 
boundary exists. 


Central South China 
Sea 
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(the region between the Dangerous Ground and the Mindoro and Balabac Straits), the 
overall performance is much better than the PSW model. 

Summer: During the summer monsoon season, both models show the surface 
currents flowing northward in response to the southwest monsoon. The PSW model 
demonstrated a good performance in the Karimata Strait and the Gulf of Thailand. 
However, it fails in the Formosa and Luzon Straits, the central South China Sea, the region 
around Palawan, and the coast of Vietnam. The PSW model results show that the waters 
from the Kuroshio Current enter the South China Sea through the Luzon Strait, and then 
flow northward into the Formosa Strait. However, the NR observations show that waters of 
the Kuroshio Current do not enter the South China Sea in this season. In the central South 
China Sea, the PSW model results show that currents flow eastward from the coast of 
Vietnam to the region around Palawan, which are oppositely directed compared with the 
NR observations. Moreover, off the coast of Vietnam, the PSW model simulated currents 
do not show jet branching flowing northeastward from the Vietnam coast to the Luzon 
Strait, which was reported by Wyrtki (1961). However, except for the region around 
Palawan, the PSCS model results match the NR observations very well. The PSCS model 
results show that the surface currents flow eastward towards the Palawan in the region 
around Palawan, which are oppositely directed compared with the NR observations. 

Winter: During the winter monsoon season, two models show that the surface 
currents flowing southward in response to the winter monsoon winds. The PSW model 
shows good performance along the coast of Vietnam and around the Karimata Strait. It fails 
in the Formosa and Luzon Straits, the central South China Sea, the region around Palawan, 
and the Gulf of Thailand. The PSW model results show that a large amount of Kuroshio 
Current flows into the South China Sea through the Luzon Strait, and then flows out of the 
South China Sea through the Formosa Strait. However, the NR observations show that a 
small amount of Kuroshio Current flows into this region through the Luzon Strait, and a 
large amount of water flows into the South China Sea through the Formosa Strait. In 


addition, in the central South China Sea, the PSW model results show that currents flow 
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westward from the western coast of Luzon and spread all over the central South China Sea, 
but the NR observations show a current boundary exists in the central South China Sea. In 
the Gulf of Thailand and the region around Palawan, the PSW model simulated surface 
currents are oppositely directed compared with the NR observations. Except for the central 
South China Sea, the PSCS model matches the NR observations very well. In the central 
South China Sea, the PSCS model results do not show the current boundary (which is 
reported by Wyrtki (1961)). 

The comparisons of current speeds between the NR observations and the PSCS 
model results are made off the coast of Hainan, off the coast of Vietnam, and the Gasper 
Strait in both summer (Table 7) and winter (Table 8). The comparisons show very good 
agreements in both summer and winter seasons. In summer, they both show that surface 
current speed is about 0.5 m/s off the coast of Hainan, more than 1 m/s off the coast of 
Vietnam, and about 0.7 m/s in the Gasper Strait. In winter, they both show about 0.5 m/s 
off the coast of Hainan, about 0.7 m/s off the coast of Vietnam, and 0.3 m/s in the Gasper 


Strait. 


TABLE 7: COMPARISON OF CURRENT SPEED BETWEEN THE NR 
OBSERVATIONS AND THE PSCS MODEL RESULTS (JUNE) 






TABLE 8: COMPARISON OF CURRENT SPEED BETWEEN THE NR 
OBSERVATIONS AND THE PSCS MODEL RESULTS (DECEMBER) 
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The PSCS model improves upon the performance of the PSW model for the 
simulation in the South China Sea. However, there are some discrepancies between the 
PSCS model results and the NR observations. The discrepancies occur in the region around 
Palawan and the central South China Sea during the summer and winter monsoon seasons, 
respectively. In summer, the NR observations show that water from the Sulu Sea enters the 
South China Sea spreading into the Dangerous Ground, but the PSCS model results show 
an eastward water flow toward Palawan. In winter, the NR observations show that a shear 
zone is present in the central South China Sea, but the PSCS model does not show the same 
pattern in this area. These might be caused by the complex bottom topography of the 


Dangerous Ground, the artificially enclosed open boundaries, and the data interpolation. 
B. TEMPERATURE AND SALINITY 


1. PSCS Model Results 


The PSCS model simulated surface temperature and salinity distributions for the 
summer and winter monsoon seasons are shown in Figures 16 and 17, respectively. We can 
clearly see that the warm water moves northward during the summer monsoon season and 
that cool water spreads southward during the winter monsoon season (Figure 16). In 
response to the change of monsoon-driven surface circulations, the water properties near 
the open boundaries also show strong seasonal variations. The East China Sea, with cool 
and high salinity waters, is an adjacent sea which can supply different water masses for the 
South China Sea. We also observe a seasonal variation of the warm, saline water in Figures 
16 and 17, northward advance during summer and southward retreat during winter. The 
surface salinity in the South China Sea is affected by the water of the adjacent seas which 
is carried by the currents and surface dilution. 

Summer: During the summer monsoon season (June) (Figure 16A), in response 
to the southwest monsoon (Figure 14A), warm water moves northeastward and there is 


upwelling of cooler (about 24°C) water along the coastlines of Vietnam and Hainan. Large 
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Figure 16A :Simulated Summer (30 June) Surface Temperature (°C) 


Distribution; Interval = 1°C 
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Figure 16B :Simulated Winter (10 December) Surface Temperature (°C) 
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Figure 17A : Simulated Summer (30 June) Surface Salinity (psu) Distribution; 


Interval = 0.25 psu 
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Figure 17B :Simulated Winter (10 December) Surface Salinity (psu) Distribution; 


Interval = 0.25 psu 


horizontal temperature gradient occur in the upwelling regions. The highest temperatures 
(about 30°C) occur in the Gulf of Thailand. Other areas of warm water are off the coast of 
Borneo and the Sunda Shelf. On the other hand, the temperature distributions are 
homogeneous (~28°C) in most of the northern parts and the other areas of the South China 
Sea. The surface temperature in the narrow area between the Asian coastline and the 28°C 
isotherm, including the Formosa Strait, the Gulf of Tonkin, and the shelf off the coast of 
Vietnam, is between 24°C and 28°C. The southwest monsoon-driven coastal upwelling, 
with surface temperatures lower than 24°C, occurs in this region. 

During the summer season, the surface salinity increases northeastward from 
32.75 psu off the Borneo coast to 34.25 psu on the western side of the Luzon Strait (Figure 
17A). Large salinity gradients occur at the coast of Borneo, the Gulf of Thailand, and the 
Gulf of Tonkin. In response to the southwest monsoon, water masses from the Java Sea 
enter the Sunda Shelf and flow northward along both the western and eastern sides of the 
South China Sea. Because of the surface river runoff into the Gulf of Thailand and around 
Borneo, the surface salinity in the Gulf of Thailand is about 32.0 psu to 33.0 psu, and is 
about 32.5 psu to 33.0 psu around Borneo. Because of river runoff, less saline water 1s 
found off the southern coast of Vietnam (surface salinity is less than 33.0 psu), and off the 
Gulf of Tonkin (surface salinity between 33.0 psu and 33.75 psu). In the northern South 
China Sea, because it connects with the cooler and more saline waters from the north, the 
surface salinity is greater than 34.0 psu. In the central South China Sea, surface salinity is 
very homogeneous, about 33.5 psu to 34.0 psu. 

Winter: During the winter monsoon season (December) (Figure 16B), in 
response to the northeast monsoon winds (Figure 14D), cool waters flow southwestward 
from the East China Sea entering the northern South China Sea through the Formosa Strait. 
The isotherms are aligned northeast-southwest increasing from the Formosa Strait to 
Singapore. Coastal upwelling disappears in this season. Because of the different water 


masses mixing strongly, the largest temperature gradients occur in the vicinity of Formosa 
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Strait with the temperature ranging from 21.5 and 22.5°C. In the Gulf of Tonkin, surface 
temperatures range between 23.0 and 25.5°C. In the Gulf of Thailand surface temperatures 
are higher than 27.5°C. Maximum temperatures of 28.5°C are found off the coast of 
Borneo. 

During the winter season, surface salinities increase from 32.50 psu in the south 
to 33.75 psu in the north (Figure 17B). The cooler high salinity waters from the East China 
Sea flow with the southwestward currents entering the northern South China Sea. The 
surface salinity in the northern South China Sea is between 33.5 psu and 34.0 psu which is 
higher than the other regions of the South China Sea. Because of the river runoff, the largest 
salinity gradients occur in the Gulf of Tonkin, the Gulf of Thailand and the constriction 
between Sumatra and Borneo. In the Gulf of Tonkin the surface salinity ranges from 31.75 
psu to 33.25 psu. On the other hand, in the Gulf of Thailand, the surface salinity ranges from 
31.5 psu to 32.75 psu. In the southern and central South China Sea, surface salinities are 
homogeneous between 32.5 psu and 33.5 psu. Between the Sumatra and Borneo, the 
surface salinity gradient is larger than the other places, and it ranges from 32.75 psu to 


33.75 psu. 


2. Comparisons 
The PSCS model results show good agreement with the Levitus temperature and 
salinity data. The comparison between PSCS model results and Levitus data is described as 


follows. 


a. Temperature 

The PSCS model simulated surface temperature (Figure 16) matches the Levitus 
temperature data (Figure 11) very well. During summer (June), the PSCS model results 
show that there is cool water off the coast of Vietnam and the southwest monsoon -induced 
upwelling occurs along the coastline of Vietnam and Hainan. The Levitus data does not 
show strong coastal upwelling, but it does show that cool water occurs off the coast of 


Vietnam. However, according to Wyrtki (1961), wind-driven upwelling is found along the 
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coastlines of Vietnam and Hainan in the southwest monsoon season. During winter 
(December), both the PSCS model results and the Levitus data show that the cold waters 


move southwestward. 


b. Salinity 

The PSCS model simulated surface salinity (Figure 17) matches the Levitus 
salinity data (Figure 12) reasonably well in most parts of the South China Sea, except at the 
coastal region of the Gulf of Thailand, the Gulf of Tonkin, and around Borneo where the 
PSCS model simulated surface salinity 1s generally greater than the Levitus salinity data 
regardless of seasons. For example, at the coastal region of the Gulf of Thailand, the PSCS 
model results show that surface salinities are about 32.25 psu in summer, and 31.5 psu in 
winter, respectively. In the same region the Levitus salinity data show that surface salinities 
are about 32.25 psu in summer, and 30.75 psu in winter, respectively. This could be caused 
by the exclusion of some physical processes in the PSCS model, such as the heavy river 


runoff. 
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V. CONCLUSIONS 


The Princeton South China Sea (PSCS) model was used to investigate the seasonal 
surface circulation in the South China Sea. An analysis of the simulated surface circulation 
shows that the South China Sea exhibits a strong seasonal variation for surface circulation 
in response to the seasonal monsoons. Also, surface temperature and surface salinity which 
are studied in this simulation show that most of the surface ocean dynamics of the South 
China Sea are strongly affected by the monsoons. 

During the summer monsoon season, southwest winds prevail and the main surface 
current flows westward from the Java Sea through the Karimata and Gasper Straits, then 
flows northeastward from Singapore along the coast of Asia toward the Luzon and Formosa 
Straits. The jet current from Vietnam to Luzon is very evident in this season. The current 
is stronger in the southern South China Sea and the current speed is about 0.5 m/s off the 
eastern coast of Hainan, more than 1 m/s off the coast of Vietnam, and about 0.75 m/s 
between Singapore and Borneo, respectively. Upwelling with a temperature of 24.0°C is 
induced and found off the coastlines of Vietnam and Hainan. Warm waters flow northward 
(warm advection) during this season. The surface temperature 1s between 28°C and 29°C 
in the South China Sea. The surface salinity is between 33.5 psu and 34.25 psu in the 
northern and central South China Sea, between 32.0 psu and 33.0 psu in the Gulf of 
Thailand, and between 32.5 psu and 33.0 psu near the coast of Borneo. 

During the winter monsoon season, northeast winds prevail and upwelling disappears. 


The main surface current flows southwestward from the Luzon and Formosa Straits along 
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the coast of Asia toward Singapore, then flows eastward through the Karimata and Gasper 
Straits entering the Java Sea. This coastal current dominates most of the South China Sea. 
Another current which flows from Borneo to Vietnam is also important during this season. 
The current speed is about 0.8 m/s off the eastern coast of Hainan and 0.7 m/s off the coast 
of Vietnam. The surface temperature is about 27.5°C in the southern South China Sea and 
is between 24.5°C and 27.5°C in the northern South China Sea. The surface salinity is 
between 33.0 psu and 34.0 psu in the northern and central South China Sea, between 31.5 
psu and 32.75 psu in the Gulf of Thailand, and between 32.25 psu and 32.75 psu near the 
coast of Borneo. 

The PSCS model simulated temperature field is much closer to the observational data 
than the salinity field, especially in the coastal regions. This is due to the exclusion of the 
river runoff process in the PSCS model. Most of the PSCS model-generated results of this 
study match the NR observations quite well and provide more detailed information than the 
PSW model. However, some model discrepancies with the NR observations still exist, 
especially the winter surface currents in the central part of the South China Sea, and the 
summer surface currents near Palawan. The discrepancies may be caused by artificially 
enclosing the four open boundaries (the Strait of Malacca, Banka Strait, Mindoro Strait, and 
Balabac Strait) or by the interpolation of sparse data sets. 

Four recommendations are provided for future studies. First, a high resolution 
database should be established including the atmospheric and the oceanic features. In this 
way, errors due to the interpolation of sparse data sets can be reduced. Second, use realistic 
open boundary conditions. Third, the model physics needs to be improved (e.g., including 


river runoff). Fourth, better measurements of the bathymetry near the Dangerous Ground 
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should be made. The discrepancies between this simulation and the observations in the 
central part of the South China Sea and the region around the Palawan (which is the region 
between the Dangerous Ground and the Mindoro and Balabac Straits), may have been 


caused by the poor bottom topographical data of the Dangerous Ground. 
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